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Abstract 

Many large scale numerical simulations of astrophysical plasmas must also 
reproduce the hydrogen ionization and the resulting emission spectrum, in some cases 
quite accurately. We describe a compact model hydrogen atom that can be readily 
incorporated into such simulations. It reproduces the recombination efficiency and line 
spectrum predicted by much larger calculations for a very broad range of densities and 
temperatures. Uncertainties in hydrogen collision data are the largest source of 
differences between our compact atom and predictions of more extensive calculations, 
and underscore the need for accurate atomic data. 



1. Introduction 



Hydrogen is the most abundant element and its physics dominates the conditions 
in a wide variety of galactic and extragalactic nebulae (Osterbrock 1989). Hydrogen's 
recombination efficiency and level populations must be computed with some precision 
since its photoionization can be the controlling heating mechanism for many clouds, it 
is the dominate opacity for many forms of light, and H is often the dominant electron 
donor (Avrett & Loeser 1988; Rees, Netzer, & Ferland 1989; Netzer 1990). The line 
spectrum must be predicted with great precision if reliable abundance determinations 
are to be achieved (Martin 1988; Shields 1990; Pagel 1992). Modern spectral synthesis 
codes (see Kallman & Mushotzky 1985; Binette et. al. 1993; Ferland et al. 1995, Netzer 
1996) must incorporate complete simulations of the H atom if they are to be valid. 

A number of physical processes affects hydrogen, including radiative and three- 
body recombination, photoionization and induced recombination, and collisional and 
radiative (spontaneous and induced) transitions between levels. Complete calculations 
of the physics of the hydrogen atom are intricate themselves (Brocklehust 1971; 
Mathews, Blumenthal & Grandi 1980, Drake & Ulrich 1980; Storey & Hummer 1995), 
and certainly cannot be incorporated into larger structure and spectral synthesis codes. 
These processes are important in the Broad Line Region (BLR) of Quasars or any dense 
gas, and so must be fully simulated. 

Storey & Hummer (1995 hereafter SH), have completed full calculations of the 
hydrogen emission spectrum and recombination efficiency over a very broad range of 
temperature and density. Their calculations are for "Case A" and "Case B" conditions, 
with the 1000 lowest n levels and all I states included. These computations are likely to 
remain definitive for some time to come, and they made their results available 
electronically, along with a convenient interpolating program. Unfortunately, it is 
possible to incorporate model hydrogen atoms with only dozens, not hundreds or 
thousands, of levels, into full spectra synthesis simulations. The problem is to define a 
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compact model atom which retains the essential physics and agrees with the definitive 
calculations. 

This paper outlines a model atom that can accurately reproduce the recombination 
efficiency and line emission with only a few dozen quantum levels, and can be 
incorporated into large-scale structure simulations. The approach goes to the correct 
asymptotic limits, and reproduces accurate results in intermediate conditions. We pay 
particular attention to obtaining accurate representations of emissivities of infrared 
lines. Surprisingly, the rates for collisions within the hydrogen atom appears to be the 
limiting factor on accuracy. Finally, we compare the idealized Case B emissivity with 
results from a realistic BLR cloud model. 

2. Physical Processes, and their simulation with a compact atom 

Our model hydrogen atom consists of independent Is, 2s, and 2p states, and an 
arbitrary number of higher excited n-levels. 

Both the low-density and LTE limits, as well as the intermediate "Case A" and 
"Case B" are well defined in the literature. Case B comes closest to simulating nature, 
and we concentrate on it here. For this approximation, all transitions involving the 
ground state are ignored, as is photoionization from excited states, collisional 
ionization and excitation from the ground and first excited states, fluorescent excitation, 
and all induced processes. These are normally included in our calculations but are 
disabled in the work presented below. 

2. 1. Photoionization, radiative recombination 

Recombination to an infinite sum of bound levels must be included if the total 
recombination coefficient, and resulting ionization balance, are to be correct. 

We have fit radiative recombination coefficients for levels up to 15 using the Milne 
relation and photoionization cross sections taken from Opacity Project-related codes, as 
described in Ferland et al. (1992). The fit has the functional form 
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a(n,T) = lO F(nT) t; 1 , 



where 



<*n + c n x + e„x + g n x + i n x 
l + b n x + d n x 2 + f n x 3 + h n x l 



4 



F(n,T e ) = 



4 



(1.) 



and x = log T e . Table 1 gives the coefficients used for the above equation, which are 
valid for the temperature range 2.8 K < T e < 10 10 K. 

For quantum levels above 15 we use an asymptotic analytical formula from Allen 
(1973): 



where % n is the ionization potential of the nth level and Ei is the first exponential 
integral (Abramowitz and Stegun 1965). 

Errors in individual recombination coefficient fits are generally 1-2% with many 
less than 1% for lower excited states at temperatures of 10,000 K. For levels 1 < n < 15 
the largest errors are 5% at only the highest temperatures. Equation 2, used for levels n 
> 16, breaks down at temperatures greater than 10 8 K, however at these high 
temperatures less that 2% of the total recombinations are to levels greater than 16, so 
this does not introduce an error in the total recombination coefficient at these high 
temperatures. 

The approach given above only provides rates for levels explicitly included in our 
model atom, and thus is an underestimate of the total recombination coefficient, a sum 
to infinity. The total recombination coefficient must be accurately computed to produce 
accurate ionization representations in the low density limit. To compensate the highest 
5 levels of the atom are assigned the remainder of the recombination coefficient: 
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where 0Cb is the total Case B rate, (the sum of recombinations to n=2 and higher levels fit 
with Equation 1 and coefficients listed in Table 1), N is the total number of levels in the 
compact model atom, and oc n is the recombination coefficient to the n th level. Although 
assigning the remaining recombination coefficient to the higher levels had the desired 
effect of including all the possible recombinations, it also resulted in the overpopulation 
of these levels. This caused a level inversion and resulted in a strong maser, and so 
was unphysical. We topped off the atom with a band of levels rather than the highest 
level to minimize this effect. The result of this topping off of the atom is that the model 
atom reproduces the ocb sum and the ionization balance at low densities to very great 
precision. 

2.2. Collisional ionization, three body recombination 

Collisional ionization and its reverse process, three body recombination, are 
important since they bring the atom into LTE at high densities. Collisional ionization 
from the ground state is taken from Arnaud & Rothenflug (1985). For levels 2 and 
higher we use rate coefficients of Vriens & Smeets (1980) who give a semi-empirical fit 
between theory and experiment for hydrogenic atoms. If the atom is in an initial 
excited state \q) with energy E q = 13.6 eV /q 2 then the rate coefficient for ionization is 

= 9.56xlQ- 6 (fcr e )- 15 exp(-e g ) 

8f 3 +4.38^ 72 +1.32e 9 ' { 

where 8 q =E q l kT e is the Boltzmann factor with kT e in eV, and the equation is good for 

the values of kT e given by Vriens & Smeets (generally for temperatures greater than 10 3 
K). For reference, the definitive SH calculations use collision data from Burgess & 
Percival (1968). Chang, Avrett, & Loeser (1991) discuss the uncertainties of the collision 
data for hydrogen. For intermediate levels this is easily a factor of two. 

Figure 1 shows the total recombination coefficient for three temperatures and a 
wide range of density. The solid line presents the results with a 50 level atom, and the 
dashed line SH. At low densities the recombination coefficient is purely radiative, 
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while three-body recombination dominates the total at high density. For low 
temperatures and intermediate densities the differences can be as great as 5%. Tests 
discussed below show that the differences in our total recombination coefficient and 
those of SH are due to the different collision data assumed, and so are a basic 
uncertainty. 

2.3. n-changing collisions 

For collisional de-excitation involving the ground, first and second excited states of 
hydrogen we use data from Callaway (1994); excitation rates are calculated by the 
principle of detailed balance. Collisions between 2s and 2p are taken from Osterbrock 
(1989). For collisional de-excitation involving levels 4 and higher we again use the 
formulae of Vriens & Smeets (1980) which is available in closed form, is recommended 
by Chang, Avrett, & Loeser (1991) , and can be extended to any transition from any 
excited state. The atom is in an initial upper level \q) and de-excites to a lower level \p) 
with energy E p = 13.6 eV / p 2 and energy difference E pq = !3.6eV (l/jo 2 - l/ q 1 ), then the 
de-excitation rate coefficient is given by 



1.6 X10- 7 (kT.r 8 P /8 



kT, +r 
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pq 
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where, g =2p 2 and g = 2g 2 are the statistical weights, and defining s=\p-q\, 
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The value 9t is the Rydberg constant and kT e is in eV, f pq is the absorption oscillator 
strength given by Johnson (1972). For comparison SH use data from Percival & 
Richards (1978). Again, differences can easily be a factor of two and are a basic 
problem (Chang, Avrett, & Loeser 1991). 

Case B does not define the populations of the n=l or n=2 levels, so collisions from 
these levels are not included in the comparisons made below. 

2.4. Radiative Cascade 

This is the most difficult process to simulate, and the only one whose treatment is 
non-standard, because the branching ratios from each level depends on both density 
and temperature. Physically, captures tend to occur to higher angular momentum, I, 
levels. As electrons cascade downward they further tend to "bunch up" at the highest 
possible / values because of the Al=l selection rule. Distant collisions with slow- 
moving protons tend to distribute the electron density population according to 
statistical weight. As a result, the actual distribution of populations among the I levels 
will have a density dependence, and the branching ratios from a particular level will 
too. 

The dashed line in Figure 2 shows the Pa to H(3 ratio from SH for three 
temperatures. The lines have a common upper level and so the intensity ratio is 
proportional to the 4-3/4-2 branching ratio. The lines are observationally important 
since both are strong and can be readily detected, and the ratio has a large dependency 
on the density and temperature of the gas. There are two simple limits; the well l- 
mixed, high density limit with Pa/H(3 = 0.277 (at all temperatures) corresponding to 
the Seaton (1959) recombination spectrum, and the low density limit with 
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Pa/H(3 = 0.339 (at 10,000 K) where I mixing collisions are negligible (corresponding to 
the calculations presented by Pengelly 1964 and Martin, 1988). 

We defined a transition probability from each excited level in terms of a constant 
total Case B lifetime (the sum of the Einstein A's for the upper level minus the ground 
state transition), and a set of branching ratios which depend on density and 
temperature as the solid line in Figure 2 shows. These branching ratios were obtained 
by fitting the SH results, for transitions from the first excited state up to n =15. Above n 
=15, the branching ratios were left at their high density values given by the Einstein 
A's. 

Branching ratios from the upper level n < 7 were calculated using two linear 
interpolations in temperature, one at log n e < 4, and one at n e = n m i n , where n m in is 
defined as the minimum density of the branching ratio. For n = 4, log n m i n = 9 (see 
figure 2). At log n e > 11 (for transitions from the upper level n = 4) the branching ratio 
is constant for all temperatures since the atom has gone over to the well /-mixed limit 
(see Figure 2). For intermediate densities between the three points, log n e =4, log n e = 
n m m, and log n e = 11, we linearly interpolate the branching ratio in density. 

Table 2a gives the coefficients of the fits used to calculate the branching ratios up 
to, and including, quantum level n =7. For log n e < 4, the branching ratio can be fit with 
the equation A + BlnT e , and at log n e = n m i n in column 3 of the table the branching ratio 
can be approximated with D + ElnT e . The value C in column 6 in the table is the high 
density branching ratio, good for log n e — nmin +2. The temperature range for all fits of 
the branching ratio is 1,000 K to 30,000 K, the valid range of SH. 

Table 2b gives the coefficients of the fits used to calculate the branching ratios for 
the atomic levels n > 8. For these levels the branching ratios were assumed to be 
simpler than Table 2a (since they are more nearly constant), such that the fits only 
interpolate in temperature for log n e < 4 (n m i n is not defined). The upper level is 
denoted in column (1) of the table and the transitions are listed horizontally to the 
right. The notation is the same as in Table 2a: at log n e < 4, the branching ratio is 
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A + B\nT e , and for log n e > 11, the branching ratio is C. For intermediate densities we 
linearly interpolate in density to calculate the branching ratio. 

Errors in the branching ratios at high densities are generally less than 1%, typically 
less than 0.5% since the branching ratios go to the high density limit for all 
temperatures. Errors at low density are larger, on the order of l%-2%. The 
intermediate density range, 4 < log n e < 9, errors are still larger, typically 5%-10%, 
because of deviation from linearity of the branching ratios for highly excited levels. In 
what follows these errors will be small compared with other uncertainties. 

3. Results 

The atom has been incorporated into the radiative equilibrium code CLOUDY 
(Ferland 1996). We have left the total number of levels of the hydrogen atom as an 
option. More levels generally produce a better agreement with the SH results, but at 
the expense of longer execution times. Tests show that lower temperatures require 
more levels, because the highest level must be well within kT e of the continuum for 
three-body recombination to achieve its full efficiency. A flexible choice in the number 
of levels was a major reason for the chosen structure of the model atom. For all of the 
results shown here a 50 level atom was used. 

Figure 3a shows a contour plot of the emissivity of H(3 (in units of 47cjv/n e n p ) for 
temperatures greater than 10 3 K and the full range of densities considered by SH. The 
H(3 emissivity varies by 5 orders of magnitude for this range in conditions. Figure 3b 
shows a ratio of this work and that of SH. For nebular temperatures (5000 K - 20,000 K) 
and all densities the differences are less than 2%. For low temperatures (< 3,000 K) and 
low densities (<10 7 cm 3 ) the general agreement is to within 6%. In this limit the 
differences are due to our use of well /-mixed Einstein A's for the higher levels (see 
section 2.4). At low temperatures captures are mainly to these levels, which are not 
actually mixed at low densities. This is a basic limitation of our approach. 
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Uncertainties in the collision data (see section 2.2) are the main reason our atom 
does not agree exactly with SH at low temperatures (<3000 K) and high densities (>10 8 
cm 3 ). We quantify the basic uncertainties by modifying the collisional rate coefficients 
from Vriens & Smeets by a factor of 2, the representative uncertainty discussed by 
Chang, Avrett, & Loeser (1991). Figure 3c shows the ratio of our predictions with and 
without this scale factor, and so shows effects of these uncertainties. Results changed by 
nearly a factor of two at low temperatures and high densities, in response. We attribute 
the differences in our results and SH at low temperature and high densities to the fact 
that we use different sources for the collision data. The results for nebular conditions 
(5000 K - 20,000 K) are affected very little, and we agree with SH very well. 

Figure 4 shows our results with SH for two Balmer lines at T e = 10 4 K relative to 
H(3. At log n e = 2.0 we calculate the Balmer decrement, Ha/Hp/Hy/HS, to be 
2.866/1.0/0.473/0.262 and at log n e = 14.0, it is 3.548/1.0/0.388/0.189. The SH results 
for the same densities are 2.86/1.0/0.468/0.259 and 3.41/1.0/0.387/0.188, respectively. 

Figure 5 shows the intensities of two observationally important infrared emission 
lines at T e = 10 4 K relative to H(3. They agree with SH to typically within 1-2%, with the 
largest errors being 5%. Tests show that this accuracy is typical for lines produced by 
levels lower than 15, the highest level we attempted to take account of variable l- 
mixing. 

The Case B approximation does not include the effects of the ground and the first 
excited states by design, so these results cannot be directly applied to dense clouds. 
Among the processes which were disabled for the comparisons above were continuum 
pumping, photoionization and induced recombination from excited states, and 
collisional excitation from the ground or first excited states. For dense clouds such as 
the BLR of a quasar, these are very important processes. 

Table 3 compares Case B with this work for 10 4 K and an electron density of 
10 11 cm 3 , typical BLR conditions Ferland et al. (1992). Column 1 lists the line labels, 
column 2 are the Case B ratios (SH) with respect to H(3, column 3 are results from our 
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compact atom with all induced processes disabled. Columns 4 and 5 show the effects 
of multiplying the collision data for levels higher than 2 by a factor of 2 and 0.5 
respectively. Clearly the results are sensitive to these uncertain numbers. Column 6 
shows the results of enabling induced photoionization and recombination and 
stimulated emission for a continuum shape and an ionization parameter (ratio of 
photon density to hydrogen density) of 0.1, given by the standard BLR model from 
Ferland et al. (1992). Column 7 has collisions from the ground and first excited states 
included in the calculation. Column 8 lists the results of allowing the Lyman lines to be 
optically thin. 

Finally, the results of a complete calculation are shown in column 9. This solves for 
the energy balance and so has a depth dependent temperature (the mean is close to 10 4 
K), it again assumes a hydrogen density of 10 11 cm 3 , has a finite column density of 10 26 
cm -2 , and corresponding H line optical depths, and a non-thermal active galactic nuclei 
continuum. In the last case most hydrogen lines are optically thick, and both 
photoionization and collisional ionization from excited levels are very important. 
Clearly the line spectrum is far from Case B. This underscores the influence of the 
ground and first excited state in high density situations. 

In summary, we find that a compact hydrogen atom can reproduce quite well the 
hydrogen emission spectrum calculated with a more extensive model atom at most 
temperatures and densities. Our results underscore the points made by Chang, Avrett, 
& Loeser (1991), indicating the need for more accurate collisional rate coefficients. 

The authors wish to thank K.T. Korista and D.A. Verner for very insightful 
discussion and questions throughout this work. We also thank the referee, G. Shields, 
for his insightful comments on our manuscript. Research in Nebular Astrophysics at 
the University of Kentucky is supported by the NSF through grant AST 93 -19034, and 
by NASA with award NAG 5-3223. 
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5. Figures 

Figure 1. Total Case B recombination coefficient defined as the radiative 
recombination plus three-body recombination minus collisional ionization compared 
with SH. This work is shown as solid lines and SH as dashed lines. Three 
temperatures are shown, 5,000 K, 10,000 K, and 20,000 K. The differences are due to 
different sources of collisional ionization rate coefficients. 

Figure 2. The SH (dashed lines) predicted Pa to H(3 ratio as a function of density 
and temperature. Three temperatures are shown, 5,000 K, 10,000 K, and 20,000 K 
moving from top to bottom. Our fit to the branching ratio is shown as a solid line. 

Figure 3. Contour plots, a) shows the log of the total H|3 emissivity (in units of 
47tjv/n e n p ) for all densities and temperatures. Solid lines are 1 dex increments and 
dashed lines are 0.2 dex steps, b) is the ratio of our predicted H(3 emissivity to SH, 10% 
increments are solid lines and 2% steps are dashed lines. Part c) of the figure is the 
same as b), but with the collisional atomic data multiplied by 2, as described in the text. 
In part c) the solid lines are 20% differences and the dashed 10%. Factors of two 
changes result for some parameters. 

Figure 4. Selected ratios of Balmer lines relative to H(3 for 10 4 K and the full range 
of density. This work is shown as solid lines and the SH results as dashed lines. 

Figure 5. Two examples of important infrared lines shown at 10 4 K and the full 
range of density. This work is shown as solid lines and the SH results as dashed lines. 
Shown are a) Br y to H(3, b) Hu a to H(3. Hu a is Humphries a the transition from level 
7 to 6 in hydrogen. 
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Fits of the b ranching ratios. 
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0.2562 


-0.0041 




5 




0.9860 


-0.0668 


0.2900 


0.4378 


-0.0145 


7 


2 


7 


-0.2313 


0.0441 


0.2360 


-0.0102 


0.0238 




3 




-0.0752 


0.0244 


0.1800 


0.0889 


0.0083 




4 




0.1399 


0.0024 


0.1636 


0.1619 


0.0000 




5 




0.3215 


-0.0135 


0.1750 


0.2684 


-0.0090 




6 




0.8889 


-0.0615 


0.2450 


0.4915 


-0.0232 



Table 2b 

Fits of the branching ratios for upper levels greater than 7 





2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


8 
A 


-0.2811 


-0.0850 


0.0523 


0.1861 


0.3673 


0.7604 
















B 


0.0465 


0.0227 


0.0080 


-0.0049 


-0.0206 


-0.0518 
















C 


0.2096 


0.1565 


0.1352 


0.1322 


0.1490 


0.2175 
















9 
A 


-0.2699 


-0.0942 


0.0149 


0.1092 


0.2100 


0.3574 


0.6726 














B 


0.0436 


0.0220 


0.0098 


0.0002 


-0.0091 


-0.0214 


-0.0452 














C 


0.1923 


0.1409 


0.1181 


0.1095 


0.1120 


0.1310 


0.1961 














10 




























A 


-0.2548 


-0.1023 


-0.0077 


0.0668 


0.1332 


0.2143 


0.3400 


0.6768 












B 


0.0406 


0.0216 


0.0108 


0.0029 


-0.0036 


-0.0108 


-0.0209 


-0.0413 












C 


0.1783 


0.1293 


0.1063 


0.0955 


0.0929 


0.0984 


0.1185 


0.1807 












11 




























A 


-0.2481 


-0.1040 


-0.0170 


0.0415 


0.0840 


0.1445 


0.2123 


0.3157 


0.5575 










B 


0.0390 


0.0209 


0.0108 


0.0043 


0.0000 


-0.0058 


-0.0115 


-0.0195 


-0.0368 










C 


0.1683 


0.1206 


0.0980 


0.0861 


0.0807 


0.0812 


0.0885 


0.1082 


0.1684 










12 




























A 


-0.2370 


-0.0988 


-0.0240 


0.0272 


0.0740 


0.1039 


0.1508 


0.2043 


0.2924 


0.5170 








B 


0.0370 


0.0197 


0.0108 


0.0051 


0.0000 


-0.0029 


-0.0072 


-0.0445 


-0.0180 


-0.0337 








C 


0.1598 


0.1138 


0.0915 


0.0791 


0.0727 


0.0685 


0.0750 


0.0802 


0.1004 


0.1586 








13 




























A 


-0.2485 


-0.1017 


-0.0311 


0.0142 


0.0660 


0.0680 


0.1112 


0.1518 


0.1983 


0.2732 


0.4787 






B 


0.0375 


0.0195 


0.0110 


0.0058 


0.0000 


0.0000 


-0.0043 


-0.0079 


-0.0115 


-0.0168 


-0.0310 






C 


0.1528 


0.1082 


0.0863 


0.0738 


0.0669 


0.0632 


0.0631 


0.0662 


0.0746 


0.0941 


0.1505 






14 




























A 


-0.2237 


-0.0995 


-0.0400 


0.0019 


0.0600 


0.0610 


0.0640 


0.1128 


0.1464 


0.1811 


0.2579 


0.4421 




B 


0.0344 


0.0188 


0.0115 


0.0065 


0.0000 


0.0000 


0.0000 


-0.0051 


-0.0079 


-0.0103 


-0.0159 


-0.0281 




C 


0.1468 


0.1037 


0.0820 


0.0697 


0.0623 


0.0581 


0.0564 


0.0572 


0.0611 


0.0698 


0.0890 


0.1437 




15 




























A 


-0.2372 


-0.0955 


-0.0370 


-0.0011 


0.0560 


0.0560 


0.0570 


0.0979 


0.1152 


0.1354 


0.18100 


0.2459 


0.4315 


B 


0.0354 


0.0183 


0.0108 


0.0065 


0.0000 


0.0000 


0.0000 


-0.0043 


-0.0058 


-0.0072 


-0.0108 


-0.0152 


-0.0260 


C 


0.1420 


0.1000 


0.0790 


0.0660 


0.0585 


0.0540 


0.0520 


0.0510 


0.0530 


0.0566 


0.0655 


0.0850 


0.1380 



Table 3 



line 


Case B a 


l b 


2 c 


3 d 


Ae 


5 f 


6§ 


BLR h 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


l-TRi 
np 


















T rv 


31 5 












3A AS 


UJ.D 


HOC 


Z.JJ 


Z.OO 


? R1 


? A1 
Z.D1 


6 71 
D./ 1 


7 11 
/ .11 


? 1Q 
Z.17 


o.zo 


riy 


n R7D 


n R7D 


n RQA 


n RAR 


D AR1 1 


n aa^ 


U.DOU 


U.oyA 


us 
MO 


n ^qa 


U.J7J 


n Ada 




U.ZDO 


U.ZDZ 


n Asn 


fl 91 1 
U.Zl 1 


niu 


U.U7J7 


U.U70J 


U.U7UD 


n 1 m 

U.1U1 


n DRQfi 

U.UJ7U 


U.UJ/ O 


niw 

U.IOO 


n 1 da 


Pa 


U.z// 


n n77 
0.1/ / 


n 177 
U.Z/V 


n 177 
U.Z// 


n 177 

U.z// 


C\ 177 
U.Z// 


n 177 

U.z// 


n inc 

O.zOo 


Pp 


0.168 


0.169 


0.175 


0.160 


0.133 


0.130 


0.185 


0.140 


Bra 


0.0654 


0.0651 


0.0679 


0.0623 


0.0515 


0.0506 


0.0719 


0.0301 


Bry 


0.0351 


0.0356 


0.0347 


0.0336 


0.0220 


0.0214 


0.0465 


0.0301 


Hua 


0.0092 


0.0093 


0.0091 


0.0088 


0.0058 


0.0056 


0.0122 


0.0022 



a Case B results from SH 

b Cloudy with physical procesess disabled as described in text. 

c Same as column 3, but with collision data times 2. 

d Same as column 3, but with collision data times 0.5. 

e induced processes enabled. 

f Collisions between 2s and 2p electrons included. 

s Case B assumption turned off. 

h Simple BLR cloud as described in the text 

1 H(3 emissivity in units of erg cm 3 sec 1 , number in paratheses is the exponent. 



